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Abstract We describe a study of period changes in 59 RR Lyrae stars, using times of 

maximum brightness from the GEOS database. The work was carried out by outstanding 
senior high school students in the University of Toronto Mentorship Program. This paper 
is written in such a way that high school or undergraduate physics and astronomy students 
could use it as a guide and template for carrying out original research, by studying period 
changes in these and other types of pulsating variable stars. 

Introduction 

The sun is a star, one of hundreds of billions in our Milky Way galaxy, which is one of 
tens of billions in the universe. Our sun makes life possible on Earth, and other stars may 
nourish life on their own planets. As stars radiate energy, however, they gradually exhaust 
their fuel, leading to change and eventual death. 

A star is an immense sphere of gas, about 3/4 hydrogen, 1/4 hehum, and 2% heavier 
elements, typically a million times larger and more massive than Earth. Throughout most 
of its life, it generates energy in its hot, dense core by thermonuclear fusion of hydrogen into 
helium, currently at the rate of 4 x 10^® W, in the case of the sun. Gradually the hydrogen 
in the core is depleted. The structure and properties of the star change. It expands into a 
red giant star. The aging of the star takes billions of years, but it has direct consequences 
that can be observed, measured, and analyzed by students! 

The Physics of Stars 

The structure of a star is governed by simple laws of physics. Deeper in the star, the 
pressure of the gas is greater, because of the greater weight of the overlying layers. Since the 
gas generally obeys the perfect gas law, the temperature and density also increase inward. 
The outward thermal gas pressure gradient balances the inward pull of gravity on the gas. 

Heat always flows from hot to less hot, so energy must flow outwards, either by radiation 
or convection, depending on the transparency of the gas to radiation. If the star is to remain 
stable, the outward-flowing energy must then be replenished by energy generation in the core; 
otherwise the core will cool, pressure will decrease, and the star will shrink. These principles, 
along with conservation of mass and energy, can be expressed as four equations. These, plus 
the input mass and composition, the gas law, and information about the transparency of the 
gas, and about the thermonuclear energy generation, can be solved on a computer to produce 
a "model" of the star - its physical properties from center to surface. The changes in these 
properties due to thermonuclear fusion of hydrogen into helium can then be incorporated 
into the calculations to produce an "evolutionary sequence of models". These predict how 
the star will age with time. They can be tested in many ways. 

Pulsating Stars 

At various stages in a star's life, its properties make it unstable to pulsation - in-and-out 
(radial) vibration being the simplest kind. The outflowing radiation interacts with the gas in 

the outer layers of the star, converting radiant energy into mechanical energy. The period P, 
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Figure 1: Light curves of an RRab (top) and an RRc (bottom) star - apparent magnitude 
(brightness increasing vertically) versus time, in units of cycles measured from maximum 
brightness. RRab stars are pulsating in the fundamental mode, and have sharp maxima. 
RRc stars are pulsating in the first overtone mode, and have rounded maxima. 



of vibration of a string (such as a guitar string) depends primarily on its length. Specifically, 
P ~ R^'^. As a result of the pulsation, the star brightens and fades from maximum to 
minimum brightness and back every P days. It is a variable star'-'^. 

A graph of the brightness of the star versus time is called a light curve (Fig. 1). It plots 
the "apparent magnitude", which is proportional to -2.51og(brightness), versus time (usually 
the "Julian Date"^) or the phase, which is the fraction of a cycle elapsed since maximum 
brightness. 

Period Changes in Pulsating Stars 

As a pulsating star ages, it initially expands slowly, so the period lengthens and the 
separation of the times of maximum brightness gradually increases. [Later in its life, the 
star may contract, in which case the period would decrease.] Specifically: (dP/dt)/P ~ 
1.5(dR/dt)/R. This process of period change, though slow, is observable because its effects 
are cumulative. Consider the analogy of a perfect clock A, and a less perfect clock B that 
runs one second more slowly every day, starting with the day when it runs at the correct 
rate. The cumulative error in clock B, after 1, 2, 3, 4 and 5 days, relative to clock A, is: 0, 
1, 1+2=3, 3+3=6, 6+4=10 seconds etc. The cumulative error increases as the square of the 
elapsed time. The possibility of watching stars slowly age by measuring the change in their 
pulsation period was recognized almost a century ago. 

Cepheids, named after the star delta Cephei"^, are an important class of pulsating stars. 
Their periods, which range from 1 to 100 days, are correlated with their average power, so 
their distance can be estimated from the inverse-square law of brightness by comparing their 
average apparent brightness with the power deduced from their pulsation periods. Polaris 
(the North Star) is the best-known star in the northern sky. It is a Ccphcid with a period 
of 3.96925... days, which is increasing by 4.5 seconds per year as the star ages. 



Figure 2: Coauthors Morenz, Meema-Coleman, and MacNeil with their poster at the UTMP 
reception and research fair, attended by mentors and mentees, family and teachers. 

Over the decades, co-author Percy's research program on variable stars and stellar evo- 
lution has been enriched by the participation of undergraduates, and of outstanding se- 
nior high school students through the University of Toronto Mentorship Program (UTMP). 
These students develop and integrate their science, math, and computing skills, motivated 
by doing real science, with real data^. Co-authors MacNeil, Meema-Coleman, and Morenz 
were participants in the UTMP in 2011. They met with Percy regularly; they produced a 
poster paper, on this project, for the end-of-year Research Fair and Reception. MacNeil and 
Meema-Coleman also developed and gave an engaging astronomy outreach program for 70 
grade six students at a local inner-city public library. 

This Project 

In 2011, for the UTMP, we decided to use pulsation period changes to study the evolution 
of RR Lyrae stars^'^ (the nomenclature of variable stars, leading to names such as RR Lyrae, 
is explained on an excellent "stars" website^ maintained by Professor Jim Kaler). They are 
old sunlike stars, with pulsation periods of about half a day, which have exhausted the 
hydrogen fuel in their core, swelled into red giants, lost about half their mass, and are now 
fusing helium into carbon. Since helium fusion is less efficient than hydrogen fusion, and 
because the star is now emitting about 100 times the power of the sun, it fuels the star for 
only a few million years, before the star again exhausts its fuel and again swells into a red 
giant. 

There are complications to be expected. Some RR Lyrae stars (designated RRab) pulsate 
in the fundamental (simplest) mode, others (designated RRc) in the first overtone, and a 
few (designated RRd) in both. Some RR Lyrae stars also show the Blazhko ejfecf, a long- 
term (weeks) variation in the amplitude and shape of the light curve which may be due to 
interference between the in-and-out pulsation mode, and a more complicated "non-radial" 
mode of pulsation. 

Data 



dally by skilled amateur astronomers. The American Association of Variable Star Observers^ 
is the world leader in coordinating such measurements; it celebrated its centennial in 2011. 
Amateur astronomers and students can make these measurements with a small telescope 
and the unaided eye, though most of the measurements are now made with CCD (charge- 
coupled device) digital cameras. We used times of maximum brightness in a database^" 
maintained by the Groupe Europeen d'Observation Stellaire (GEOS), including times from 
previous publications, from the AAVSO and other such organizations, and from robotic 
telescopes^^. The database includes about 50,000 times of maximum of over 3,000 stars, 
covering up to a century. There were 59 stars with data that were long and continuous 
enough for straight-forward analysis. 

Analysis 

For any star, the observed times of maximum brightness t(i) are compared with the 
predicted times t(0) -|- NP where t(0) is an initial time of maximum, P is the assumed 
period in days, and N is an integer - the number of cycles between t(0) and t(i). t(i) - (t(0) 
+ NP) is called (0-C), observed minus calculated. A graph of (0-C) versus t - an (0-C) 
diagram or OCD - will be a straight line if the period P is constant; its slope will be positive, 
zero, or negative, depending on whether P is greater than, equal to, or less than the true 
period. This slope can therefore be used to correct or refine the existing period of the star. 
If P is changing hnearly (i.e. P = P(0) -|- kt where k is a constant), then the OCD will be a 
parabola, opening upward or downward, depending on whether k is positive or negative. A 
parabola can be fitted to the OCD using least-squares software, such as the fitting function 
in spreadsheet software such as Excel. 

A complication occurs if there are long gaps between the t(i): the number of cycles in 
the gap is uncertain. If the gaps are not too long, and if the measurements before and after 
the gaps are reasonably continuous, it should be possible to infer the number of cycles in 
the gap by extrapolation of the existing measurements. 

The database was scanned for stars having a sufficient number and continuity of t(i). 
Periods and t(i)s were entered into a spreadsheet; values of C and of (0-C) were determined; 
OCDs were plotted. If their slopes were non-zero, the existing periods were modified. 

We initially classified the OCDs as approximately (i) parabolic (P); (ii) linear (L); (iii) 
wavy (W); (iv) otherwise peculiar (?) - not parabolic or linear; or (v) too scattered (S) to 
classify. We later decided that parabolic-linear was a continuum in the sense that "apparently 
linear" corresponded to a parabola with very small curvature, and that "wavy" was a special 
case of "pecuhar" . 

The accuracy of the t(i) is greater for the RRab stars, because they have sharp-peaked 
light curves, whereas the RRc have more sinusoidal ones (Fig. 1). The accuracy is generally 
also greater if the t(i) has been measured with a CCD camera, rather than with the unaided 
eye or photography. The accuracies are typically a few minutes. 

A characteristic time scale, for period change and therefore for change in the size of the 
star, is the period divided by the rate of change of the period P/(dP/dt), in appropriate 
units, i.e. it is the time required for the period to change by P. We designate the time scale 
by T (tau). This is approximately the length of time that the RR Lyrae phase of evolution 
would last, though that time more correctly depends on how long the star's properties cause 
it to be unstable to pulsation. Alternatively, the time scale could be defined as the time 
required for the period to change by a factor of 2, or e. These time scales apply to the few 
decades which are covered by our data, though evolutionary models predict that the rate of 
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Results 

Tables 1 and 2 list the stars analyzed: their name, period (day), (0-C) diagram type, 
RR Lyrae subtype, [Fe/H], dP/dt in seconds per 100 years, and the characteristic evolution 
time r = P/dP/dt in millions of years. [Fe/H] is the logarithm of the ratio of a star's iron 
abundance to its hydrogen abundance, compared with that of the sun; it is used as a measure 
of the abundance of elements heavier than helium. 

1. Of the 59 stars that were analyzed, 10 were too scattered (S) to classify, 27 were 
parabolic (P), 14 were linear (L), P + L = 41 in total; 3 were wavy (W) and 5 were 
otherwise peculiar (?), W + ? = 8 in total. 

2. The two RRd stars were both among the stars with too much scatter for analysis; their 
period changes must be studied in more sophisticated ways. 

3. The Blazhko-effect stars showed no particular tendency for scatter, though 5 of the 18 
had peculiar OCDs. 

4. Of the RRab stars whose (0-C) diagrams could be fit with parabolas, 15 had increasing 
periods (median r of 10 million years) and 14 had decreasing periods (median r of 6.7 
million years) i.e. there were approximately the same number increasing as decreasing, 
and the rates of evolution were similar. 

5. The RRc stars have sinusoidal light curves. The times of maximum are more difficult to 
measure than for the RRab stars, whose light curves have sharp maxima, so they were 
not studied by Le Borgne et al.^^. Of the 14 RRc stars in our sample, 4 had scattered 
OCDs, 2 had pecuhar OCDs, 5 showed increasing periods (median r of about 20 miUion 
years), and 3 showed decreasing periods with a large scatter in the rate. The median 
r suggests rather slow evolution, compared with the RRab stars. The sample size, 
however, is small. 

Two examples of OCDs are shown in Figs. 3 and 4. RR Leo has a parabolic OCD, 
corresponding to a period increase of 2.7 seconds per century, and RR Gem has a peculiar 
one, neither parabolic nor linear. A few OCDs seem to show abrupt pulsation period changes, 
which are difficult to explain theoretically. Various explanations have been proposed for these 
peculiar OCDs, including "burps" of convection^^, deep in the star, or magnetic cycles in 
their outer layers. 

Discussion 

It is not easy to compare our results in detail with the predictions of evolutionary models^^ 
because the sign and magnitude of the predicted changes depends strongly on the exact mass 
and composition of the star, especially the abundance of the elements heavier than helium 
(denoted Z). The value of [Fe/H] (see Tables 1 and 2) is used as a proxy for Z. Indeed, the 
stars in our sample have a wide range of values of [Fe/H], and probably of masses (though 
most are probably between 0.4 and 0.6 times the sun's). The magnitudes of the observed 
period changes, and the values of r are reasonably consistent with the models, though the 
models predict many more period increases than decreases. Sec references 11-13 for more 
detailed comparisons between observations of period changes, and predictions of evolutionary 
models. 

Pulsation period changes have been studied in other kinds of variable stars, with a variety 



0-C vs. Julian Date; RR Leo 
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Figure 3: The parabolic OCD for RR Leo, corresponding to a linear period increase of 2.7 
seconds per century. 



pulsate with periods of several hundred days. In these stars, the pulsation drives off the outer 
layers of the star into space, leaving a white dwarf corpse, a million times denser than water. 
The pulsation of these stars is complicated by random cycle-to-cycle period fluctuations 
of a few percent, leading to wave-like (0-C) diagrams. Times of maximum brightness are 
available for analysis on the AAVSO website^^. 

Eclipsing binary stars are pairs of stars in mutual orbit in which one or both stars 
periodically eclipses the other star, resulting in a brief decrease in brightness. (0-C) analysis 
of the times of minimum brightness, also available on-line^^ provides information on mass 
transfer or loss in the system. 

The period changes in Cepheid pulsating stars are particularly important for testing 
evolutionary models, and have been extensively studied by Leonid Berdnikov, David Turner, 
and their colleagues^^. High school (and university) students can observe the variability 
of Cepheids using the unaided eye^^ or a camera^^, and verify the slow period changes in 
these stars. Students' observations of these or other variable stars can be entered into the 
AAVSO International Database, which now contains over 20 million observations, for other 
researchers to use. 

The detectability of linear period changes, due to the aging of the star, increases as 
the square of the timespan of the data, so it is worthwhile for observers to continue their 
measurements, and for analysts to repeat the (0-C) analysis as longer datasets become 
available. For instance: the period changes in some of the RRab stars in our sample were 
studied a few years ago^^ but the analysis should be repeated periodically. We believe that 
this type of analysis provides an engaging education experience for undergraduates or for 
senior high school students; the physical and mathematical concepts are simple, data is 
available on-line and, with it, students can actually observe the effects of the slow aging of 
the stars. 



Conclusions 
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0-C vs. Julian Date; RR Gem 
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Figure 4: The non-parabolic, non-linear OCD for RR Gem. The OCD is not explainable as 
due solely to evolution, but must be due to one or more other processes. 

observe and measure the evolutionary changes in the stars, even though the characteristic 
times, for the changes, are millions to billions of years. In this way, they can make meaningful 
use of their science and math skills, motivated by doing real science, with real data. 
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RR Lyr 





56683295 
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-1 


14 
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37733834 
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RU GVn 
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-1 
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RZ Gep 
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-1 
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RZ Lyr 





51124591 
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ab 


-1 


69 


-1.578 


-2.71 


S Gom 





58658810 


p 


ab 


-1 


91 


-0.892 


-5.50 
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36733833 


p 


ab 


-0 


24 


0.185 


16.6 


SS Tau 





36991119 


L 


ab 






-1.270 


-2.44 


ST Boo 





62229368 


L 


ab 


-1 


72 


0.348 


15.0 


ST Leo 





47798421 


L 


ab 


-1 


17 


-0.416 


-9.61 


SU Dra 





66042067 


P 


ab 


-1 


39 


0.534 


10.3 


SW And 





44226806 


P 


ab 


-0 


24 


-1.446 


-2.56 


SW Aqr 





45930290 


L 


ab 


-1 


63 


0.067 


57.4 



Table 1: Rates of period change, and characteristic evolution times for RR Lyrae stars. 
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56967136 


P 


ab 


-1 


12 


0.353 


13.5 


SX Aqr 





53571243 


P 


ab 


-1 


87 


-0.686 


-6.53 
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43605526 


P: 


ab 




— 


-0.894 


-4.08 


SZ Hya 
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? 


d 




— 


— 


— 


TU UMa 





55765749 


L: 


ab 


-1 


05 


-0.160 


-29.1 


SX UMa 
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? 


c 


-1 


82 


— 


— 


AA Aql 





36178731 


P: 


ab 


-0 


20 


0.032 


94.6 


AH Cam 
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B 


ab 




— 


-4.100 


-0.75 


AQ Cnc 





54851950 


P: 
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— 


-0.353 


-13.0 


AR Per 





42554959 


P 


ab 


-0 


17 


0.250 


14.2 


AV Peg 
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P 


ab 





02 


1.537 


2.13 
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47632365 


P 


ab 


-0 


15 


0.499 


7.99 


AC And 
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S 


d 


-1 


16 


— 


— 


AN Ser 
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s 


ab 


-0 


07 


— 


— 


AQ Lyr 
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S: 


ab 
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?: 


ab 


-1 
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S 


ab 


-1 
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57514160 


S 


ab 


-0 


85 






BC Dra 
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L: 


ab 


-2 


00 


-0.777 
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BD Dra 
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P: 


ab 
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1.16 


TV Boo 
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P 


c 






0.087 


30.1 
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L: 


c 






0.248 


10.6 


ST CVn 





32904500 


W 


c 
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S 


c 










VZ Dra 





32102829 


P 


c 






-1.441 


-1.93 


LS Her 





23080777 


P: 


c 






0.052 


42.4 


TV Lyn 





24065125 


P: 


c 






0.091 


22.9 


DH Peg 





25551056 


L: 


c 






-0.037 


-60.3 


SS Pse 





28779248 


P 


c 






-0.207 


-10.2 



Table 2: Rates of period change, and characteristic evolution times for RR Lyrae stars 
(continued) . 



